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Abstract. It has been proposed that the observed dark energy can be explained 
away by the effect of large-scale nonlinear inhomogeneities. In the present paper we 
discuss how observations constrain cosmological models featuring large voids. We start 
by considering Copernican models, in which the observer is not occupying a special 
position and homogeneity is preserved on a very large scale. We show how these models, 
at least in their current realizations, are constrained to give small, but perhaps not 
negligible in certain contexts, corrections to the cosmological observables. We then 
examine non- Copernican models, in which the observer is close to the center of a very 
large void. These models can give large corrections to the observables which mimic an 
accelerated FLRW model. We carefully discuss the main observables and tests able to 
exclude them. 

PACS numbers: 95.36.-t-x, 98.62.Sb, 98.65.Dx, 98.80.Es 
1. Introduction 

A striking feature of the late universe is the abundance of large voids seen both in 
galaxy redshift surveys [1, 2, 3, 4, 5] and large-scale simulations [6]. Voids dominate 
the total volume of the universe while matter is mostly distributed into a filamentary 
structure dubbed as the cosmic web. Photons, therefore, travel mostly through voids 
and it has been asked by various authors if large-scale nonlinear inhomogeneities can 
sizeably alter the observables as compared to the corresponding homogeneous and 
isotropic Friedmann-Lemaitre-Robertson- Walker (FLRW) model. The study of the 
differences between the "optical universe" and the homogeneous counterpart started 
perhaps with Zel'Dovich [7] in the 1964 and continued in the following years by, for 
example, Refs. [8, 9, 10, 11, 12, 13, 14]. 

After 1998, when the dark energy became observationally necessary [15, 16], it was 
then asked (see, for example, [17, 18, 19, 20, 21] and references therein) if the corrections 
could be large enough so to have a paradigm shift in which the dark source is no longer 
needed. The large-scale structures became indeed nonlinear recently, exactly when a 
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primary dark energy is supposed to start dominating the energy content of the universe 
and cause acceleration. In this paper we will review some of the contributions to this 
topic and how they are constrained by experimental data. The common thread will be 
the modeling of voids. 

The effect of large-scale inhomogeneities on cosmological observables has been 
studied in the literature following two non- exclusive approaches. With the first it is 
assumed that one can compute observational quantities as in an effective homogeneous 
and isotropic background metric, whose evolution equations however are modified, as 
compared to the usual Friedmann equations, by the presence of large-scale structures 
due to the nonlinear nature of gravity; this topic will be covered in detail by other 
reviews of this Special Issue [22, 23, 24, 25, 26, 27, 28, 29]. The second approach, on the 
other hand, focuses on the observational properties of the universe: even in the case in 
which one ignores the nonlinear effects of gravity on the overall background evolution, 
there are still potential differences due to the fact that photons in a inhomogeneous 
universe propagate differently than in an FLRW model. The former effect has been 
called sometimes "strong" backreaction and the latter "weak" backreaction (see Ref. [30] 
for more precise definitions). As will be clear later on, the models discussed in this paper 
address the weak backreaction issue, as strong backreaction is by construction small in 
these examples. 

This paper is organized as follows. In Section 2 we start by considering Copernican 
models, in which the observer is not occupying a special position and homogeneity is 
preserved on a very large scale. We then examine in Section 3 non-Copernican models in 
which the observer is close to the center of a very large void. We give our conclusions and 
discuss future prospects in Section 4. In Appendix A we briefly introduce the spherically- 
symmetric LTB metric and in Appendix B we argue that strong backreaction effects are 
small for realistic exact Swiss-cheese models. 

2. Copernican Models 

Within Copernican models the observer is not occupying a special position and the 
corrections to the observables come from the cumulative effect of many inhomogeneities. 
In this Section we will consider two particular configurations, namely, the so-called 
"swiss- cheese" and "meatball" models. 

2.1. Swiss- Cheese Models 

The first swiss-cheese models go back to Einstein and Strauss and consist of 
Schwarzschild holes embedded into an Einstein-de Sitter (EdS) background [31], and 
the idea was to hide matter from the observer by confining it into the Schwarzschild 
masses at the center of the swiss-cheese holes [11]. This idea will be developed with the 
meatball models in the next Section, while here we will focus on the modelling of the 
voids by considering a swiss cheese with LTB holes. Lemaitre-Tolman-Bondi models 



Observational constraints on inhomogeneous models without dark energy 




Figure 1. Sketch of the Swiss-cheese model of Ref. [38, 39]. The shading mimics the 
density profile: darker shading implies larger densities. The uniform gray is the EdS 
cheese. The photons pass through the holes as shown by the arrows and are revealed 
by the observer placed in the cheese. 



[32, 33, 34] are spherically symmetric dust solutions of Einstein's equations with pure 
radial motion and without shell crossing (see Appendix A for a brief introduction). 

The solutions we are considering involve only pressureless matter. However, this 
should not affect the modelling of the voids at late times as pressure is negligible within 
them (see [35, 36] for models which contain also radiation). Moreover, the spherical 
approximation should not be too unrealistic as voids tend to evolve towards a spherical 
configuration. It can indeed be shown with the top-hat void model that the smaller 
axis of an underdense ellipsoid grows faster as compared to the longer ones, with the 
consequence that voids become increasingly spherical as they evolve [37]. The exact 
opposite happens to overdense ellipsoids which tend to form filamentary and pancake- 
like structures. On the other hand, simple spherical voids neglect the possible presence 
of substructure and the interactions between the underdensities. 

The general setup of an LTB patch is of an inner underdensity ( "hole" ) surrounded 
by a compensating overdensity ("crust"), necessary to match the EdS background 
("cheese"): see, for example, the sketch of Fig. 1. A physical picture is that, given 
an EdS sphere, all the matter in the inner region is pushed to the border of the sphere 
while the quantity of matter inside the sphere does not change (note that this is valid 
only for holes matched at a finite radius). More precisely, in order to match exactly 
the two metrics the average density of the inhomogeneous region has to be equal to the 
density of the external FLRW metric, f With the density chosen in this way, an observer 
outside the hole will not feel its presence as far as local physics is concerned (this does 
not apply to global quantities, such as the distance-redshift relation for example). So 

| A subtlety is that the average density has to be defined with a volume element which does not include 
the curvature term (see the discussion about the effective mass function F(r) in Appendix A and also 
Ref. [34, 40] ) . 
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the cheese is evolving as an EdS universe while the holes evolve differently. In this way 
we can imagine putting in the cheese as many holes as we want, even with different 
sizes and density profiles, and still have an exact solution of Einstein's equations (as 
long as there is no superposition among the holes and the correct matching is achieved). 
This property is due to spherical symmetry and it is qualitatively analogous to Gauss' 
theorem in electrodynamics or to Birkhoff 's theorem for black holes. While this allows 
to solve the equations exactly, it eliminates on the other hand any important "strong" 
backreaction [30] effect in the FLRW region, which keeps evolving exactly as in a pure 
FLRW model no matter how the hole is modeled. For more details see Appendix B 
where we argue that strong backreaction effects are small for realistic exact Swiss-cheese 
models (see also [40] for a discussion about strong backreaction in LTB models which 
converge asymptotically to FLRW). For this reason it may be of crucial importance to 
study metrics which do not have the property of being exactly glued to FLRW, while 
still having homogeneity and isotropy on large scales. So far, to our knowledge, all the 
models in the literature are based on LTB or its generalizations, and thus still preserve 
this property, even in absence of spherical symmetry (see, for example, Refs. [41, 42, 43] 
for models based on the Szekeres metric; see Ref. [44] for models which only assume 
statistical homogeneity and isotropy; see Ref. [45] of this Special Issue for an extensive 
discussion of inhomogeneous models relevant for cosmology). 

LTB swiss-cheese models have been studied in the literature with different 
configurations and sizes, analytically and numerically. In order to interpret cosmological 
observations the key quantities to be computed for a given source are redshift and 
distance, as for instance in any supernova analysis. Inhomogeneities affect both 
observables through redshift and lensing effects, which we now discuss. 

2.1.1. Redshift Effects. When a photon crosses a hole, it gains an additional redshift 
5z with respect to the usual FLRW redshift. It is, however, natural to expect a 
compensation on the redshift effects between the ingoing and outgoing geodesic paths 
due to the spherical symmetry. Moreover, because the LTB metric is matched to the 
background metric, it can be shown that there is a compensation already on the scale 
of half a hole. The effect on a photon has been studied exactly in the full nonlinear 
LTB metric [38, 46, 47, 48] both numerically and using analytical approximations. It 
is, however, possible to capture the qualitative meaning of the various leading effects by 
means of a perturbative approach which should be valid if the size of the hole is much 
smaller than the FLRW horizon and could perhaps give a better idea to the reader of 
the physics behind the exact results. Of course this discussion does not substitute the 
exact nonlinear one, which is anyway necessary when the LTB patch extends to sizes 
comparable to the Hubble horizon. 

The correction to the redshift can be understood via the well-known approximate 
expression [49] for the redshift (or equivalently for the temperature of photons if they are 
blackbody distributed) in a given direction e*, which is expressed here in the Newtonian 
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gauge and is valid for small velocities and potentials: 

5z 5T f E <9$ 

$E-&o-v % j S e i + v % ei-2 dr — . (I) 



1 + z T — w , -j^ dr 

In Eq. (1) the first two terms give the so-called Sachs- Wolfe effect, due to the difference 
in gravitational potentials $ at Emitter and Observer, the third and fourth are the 
Doppler effects due to the velocity v 1 of the Emitter and the Observer and the last term 
is the Integrated Sachs Wolfe (ISW) effect, due to the change in (conformal) time r of 
the potential along the line of sight. Each of these terms can be explicitly computed 
in the LTB metric in the "small-u" approximation, which has been extensively studied 
in [46, 47]. 

As long as the proper radius of the hole Iq is much smaller than the size of the 
FLRW horizon Zhor, we can easily understand which are the dominant terms. If the 
observer or the source occupy a generic position, the dominant terms are the velocities, 
which can be shown to be of order /o/^hor- If we put the source at the boundary and the 
observer at the centre (or viceversa), the velocities vanish, because of the exact matching 
to FLRW, and there are only two nonzero terms: $o(b) and the ISW effect, of which the 
dominant one is $o(e) which can be shown to be of order (ZoAhor) 2 - Finally, if we put 
both the observer and the source on the boundary, the only non-vanishing contribution 
is the ISW, which is further suppressed and goes as (ZoAhor) 3 - This last effect is also 
known as the Rees-Sciama effect [50] and comes from the fact that the hole and its 
structures evolve while the photon is passing, and it has been confirmed analytically in 
LTB models by Ref. [47, 48], in which more details can be found. 

The suppression of redshift effects due to the matching of the LTB metric to the 
EdS background can be also qualitatively understood by noting that dz/dr = H oc p = 
Psds + dp- Because the density profile satisfies (5p) = 0, in its journey from the center to 
the border of the hole (or viceversa) the photon will see a (H) ~ H EdS and there will be 
compensation for z' . This can also be seen by line averaging the longitudinal expansion 
rate Hl (see Appendix A for definitions) from the center to the boundary of the hole 
[38]: 

j;°drH L Y' _ Y 



(fh) ^ ■ - - , (2) 

r=r 

where Y is a r-dependent generalization of the scale factor for LTB metrics, ro is the 
comoving radius of the LTB patch (Y(r Q , t) = a(t) r = l ) and the approximation comes 
from neglecting the small curvature E and the time evolution of the density profile. The 
latter can be qualitatively interpreted as neglecting the ISW effect and the former as 
neglecting the Sachs- Wolfe effect ($ oc E, see [47]). 

In conclusion, redshift effects are large, and actually the dominant ones if one 
considers sources or observers within the holes, and especially at low redshift. This has 
been discussed in detail for the Swiss-cheese model of Ref. [47] (left panel of Fig. 2), 
and also for the so-called Onion model of Ref. [46] (see also [51]), where the universe is 
described by spherical concentric shells with positive and negative contrast (right panel 
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Log 10 (L/Mpc) r _ r0 

Figure 2. Left panel: net correction to the redshift (5z) for the Swiss-cheese model 
of Ref. [47] as a function of the size L of a hole with S s» 0.2, for an observer sitting 
outside the hole. The triangles are the numerical results, while the solid line shows 
the predicted cubic dependence (5z oc L 3 ). Right panel: redshift along the geodesic of 
a photon revealed by an observer far from the center of the Onion model of Ref. [46] . 
The solid line is the numerical solution, the dashed line is the analytical approximation 
and the dotted line is the FLRW result. 



of Fig. 2). Similar results can be seen also from Fig. 3 (left panel), which is relative to 
the Swiss-cheese model of Refs. [38, 39]. However these are mainly due a Doppler effect 
from peculiar velocities, and this should average out when considering many sources, 
except in the case of non-Copernican models, which will be addressed in Section 3 of 
this review. If we indeed consider light going through the Swiss-cheese holes and both 
observer and emitter in the cheese, then the effect from a single hole is suppressed as 
{h/Worf ■ If now light goes through a series of many voids, the effect basically sums up: 
in the extreme case of filling the universe with LTB patches, a photon may meet I^/Iq 
number of holes and this gives an overall correction of order (ZoAhor) 2 ; which would not 
be enough to explain away dark energy, for any realistic configuration with voids smaller 
than the horizon (see [52] for an extensive numerical analysis). 

2.1.2. Lensing Effects. In addition to the effects on the redshift of a source, it is 
necessary to study the effects on its distance, which can be defined equivalently as 
the angular distance d 2 A = dA/dVL or the luminosity distance d^ = (1 + z) 2 d A , 
where dA and dVt are the area of an object and the solid angle under which it is 
seen by an observer, respectively. Rather than the distances above, we will use 
the distance modulus which is defined as m(z) = 51og 10 dz,(2;)/(10pc). Large-scale 
inhomogeneities affect distance measurements because their lensing effects change the 
ratio dA/dVt. This has been computed by several authors in the exact LTB model by 
solving explicitly for the geodesies in the full LTB metric, both numerically and with 
analytical approximations [38, 39, 46, 47, 52, 53, 54]. Lensing effects can be also simply 
understood in the weak-lensing limit of small magnifications, which we briefly introduce 
as follows. For more details see, for example, Ref. [55]. 

In the weak-lensing theory the net magnification /i produced by a localised density 
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Figure 3. Left panel: change in redshift with respect to the EdS model for a photon 
that travels from one side of the five-hole chain of Fig. 1 to the other where the observer 
will detect it at present time. The vertical lines mark the edges of the holes. The plots 
are with respect to the coordinate radius r. Notice also that along the inner voids 
the redshift is increasing faster because of the higher expansion rate (z'(r) = H(z)). 
Right panel: the luminosity distance as a function of redshift for the observer of Fig. 1, 
together with the EdS and ACDM curves. Rather than Hodi(z), we show the usual 
difference in the distance modulus compared to the empty model. 



perturbation is: 



1 



^ = (l-«)2-| 7 ,|2 ^ (!-«)" > ( 3 ) 



where re is the lens convergence and we have neglected the second-order contribution of 
the shear j 8 . The convergence is due to the local matter density (see Eq. (5) below) 
and its physical effect is to magnify the image by increasing its size; the shear instead 
is due to the tidal gravitational field and its effect is to stretch the image tangentially 
deforming a circular shape into an elliptic one [55]. The shift in the distance modulus 
caused by \x then becomes: 

Am = -2.5 log 10 n ~ 5 log 10 (l - re) . (4) 

The lens convergence re can be computed from the following integral along an 
unperturbed light geodesic: 

k{z s )= ^ dr r{rs ~ r \ 2 ^, (5) 

JO r s 

where r labels the comoving distance and z s is the redshift of a light source whose 
comoving position is r s = r(z s ). The term V 2 $ is the Laplacian of the Newtonian 
potential, given by 

V 2 $ = AnGa 2 5p M = | ^ 8 M , (6) 

where a(t) is the EdS scale factor, 5m = &Pm/ Pm is the matter contrast and pu is the 
matter density. From Eq. (5) follows then 

k(z ) - -H 2 r dr ^ - r) Mr ' t(r)) m 
k{z s )- 2 H J o dr ^ . (7) 
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The previous equations show that for a lower-than-EdS column density the light is 
demagnified (e.g., empty beam 5m = — 1), while in the opposite case it is magnified. The 
weak-lens approximation is valid when the universe can be described by a Newtonian- 
perturbed FLRW metric and when k C 1, i.e., when the average matter contrast along 
the line of sight is 5 m < 1- 

As shown in Fig. 3 (right panel), Refs. [38, 39] found a large effect for the setup 
where an observer in the cheese looks through a chain of aligned voids (see Fig. 1). 
The idea was indeed that the photons we observe have likely travelled more through 
(large) voids than through dense structures. The effect builds up and it becomes large 
at high redshift z > 1. The particular configuration chosen in [38, 39] had holes of a 
radius of 350 Mpc; however, it can be shown [52] that the effect does not depend on the 
size of the voids as long as the average matter contrast along the line of sight remains 
approximately constant (see Eq. (7)). For the setup of Fig. 1, a chain of 50 holes which 
are 10 times smaller will indeed yield a lightcone column density approximately equal 
as compared to the configuration with 5 larger holes. For constant-time slices this is 
actually exactly true, since the LTB dynamics is invariant with respect to the hole 
radius once the density profile is properly scaled [38]. 

Such a result, however, can be valid only for particular directions in the sky, such 
as the one where all the centers of the voids are aligned, and not for a generic direction. 
Indeed, since the distance is a measure of the total luminosity received and the number 
of photons is conserved, the average of and therefore over the full sky must 
be identical§ to the unperturbed sky, at least in the case of weak lensing where no 
photons are lost along the path||. In other words, while some areas of the sky are 
demagnified other areas are magnified, yielding an exact compensation for the full-sky 
average, showing the "benevolent" nature of weak lensing. So, the only overall net effect 
for the full-sky average is in principle the change in the photon redshift discussed in the 
previous Section. This has been confirmed by Refs. [52, 54, 59] (see also [53]) where, 
by averaging over the lines of sight or by randomizing the positions of the holes, the 
luminosity distance has been found to converge to the EdS prediction. An additional 
relevant result is that the Rees-Sciama effect on the CMB constrains the holes of a 
swiss-cheese model to have a radius smaller than about 35 Mpc [52], thus making it 
unlikely that a photon always travelled through the central regions of the voids as in 
Fig. 1. 

It is interesting to note that, even if the average is preserved, the lensing probability 
distribution function (PDF) is modified by the presence of structures. The Swiss-cheese 
model, however, is not suitable to study the statistical properties of light propagation 
and the basic reason is that a photon always hits the surrounding overdensity before and 
after passing through the underdensity. This geometrical feature imposed by the need to 
match the EdS metric constrains the model to have a nearly gaussian PDF [52, 53, 59], 

§ This should be approximately true [13] if strong lensing events leading to secondary images and 
caustics [14] do not play a significant role as far as the full-sky average is concerned [56, 57, 58]. 
|| This can be seen by the fact that it is {8m) = and so Eq. (7) implies (k) = 0. 
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especially if the voids, as discussed above, have to be smaller than about 35 Mpc. In 
order to properly study the statistical properties of the luminosity distance in a universe 
dominated by voids, i.e., the lensing PDF, it is necessary therefore to focus on a model 
that allows photons to miss overdensities: this will be the subject of the next Section. 

2.2. Meatball Models 

The meatball model (see, for example, [60, 61] and also the lattice model of [62]) 
describes the universe as made of a collection of possibly different spherical halos^f and 
incorporates quantitatively the crucial feature that photons can travel through voids, 
miss the localized overdensities and experience a low matter column density. This is 
due to the fact that the underdensities occupy more volume than the overdensities and 
causes the lensing PDF to be skewed with a mode at demagnified values. This is of 
potential importance for the interpretation of the supernova observations, which are 
still probing much smaller angular scales than the scale at which the homogeneity is 
recovered. Moreover, mechanisms able to obscure lines of sight are generally related 
to mass concentrations, with the consequence that selection effects give a neat bias at 
demagnified values: in other words, they systematically hide matter from observations. 
Selection effects connect indeed the meatball model to the empty or partially-filled beam 
formalism started by Zel'Dovich [7] and developed later by Kantowski, Dyer, Roeder 
and other authors [8, 9, 10, 11, 12, 13, 14, 64, 65]. We will discuss here an illustrative 
example where the survival probability is modeled with a simple step-function in the 
impact parameter. 

It could be possible that lensing effects due to a skewed PDF, possibly strengthened 
by selection effects, could explain away part of the necessary dark energy, also because 
lensing effects are stronger on a EdS background as compared to the concordance model. 
To illustrate this point we compute the lensing PDF relative to the EdS model with 
h = 0.5 (as always H = 100 h km s -1 Mpc -1 ) using the turboGL code [66] based on 
the sGL method [67, 68, 69]. + In this example the meatball model has two families of 
spherical halos, each one comprising half of the total matter density. The halos of the 
first family have a mass of 10 14 /i~ 1 M and SIS (singular isothermal sphere) profile with 
radius of R = 1 h~ x Mpc, while the halos of the second family a mass of 10 17 h~ l M & and 
Gaussian profile with radius of R = 10 h^ 1 Mpc. Moreover, we toy model the selection 
effects by obscuring light that hits the halos with impact parameter smaller than 20% of 
the radius. We show in Fig. 4 the lensing PDF without selection effects (taller histogram) 
and the PDF with selection effects (shorter histogram). While the former has correctly 
a mean at zero (background) value, the latter has a mean at the (demagnified) value 
of Am ~ 0.09 mag. Moreover, the mode is in both cases at large demagnified values, 

% See Ref. [63] for a reanalysis of the SNe data which includes the lensing effects of the meatball model. 
+ The sGL method for computing the lensing convergence is based on generating stochastic 
configurations of halos and filaments along the line of sight. The turboGL code is a numerical 
implementations of the sGL method and allows to compute the full lensing PDF in a few seconds. 
See, for example, [56, 57, 58, 70] for lensing PDFs obtained through ray-tracing techniques. 
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Figure 4. Lensing PDFs in magnitudes for the meatball model described in Section 
2.2 for a source at z = 1.5. The taller (green) histogram is without selection effects, 
while the shorter (orange) histogram is with selection effects. While the former has a 
mean at zero (background) value, the latter has a mean at the value of Am ~ 0.09 
mag. The y-axis follows from the normalization J dAmP(Am) = 1. 

Am > 0.2 mag. This toy example shows, therefore, that the skewness of the lensing PDF 
together with selection effects may mimic the demagnifying effects of the cosmological 
constant. It is also interesting to note that, once selection effects have been properly 
modeled, experimental data (e.g. SNe datasets) can constrain the lensing PDF and so 
the spectrum of matter inhomogeneities [68] . 

Lensing effects alone, however, cannot explain away dark energy and the basic 
reason is that they are non- negligible only for z > 0.5. Nonetheless their effect could 
add up to other effects. They could, for example, enhance non-Copernican models (see 
next Section) or make the necessary void radius smaller [61]. Moreover, in the above 
toy model redshifts are related to comoving distances through the background metric, 
thus neglecting any redshift effects, which may be significant. They may lead indeed 
to large corrections to the observables, especially in presence of selection effects which 
may change the number of blueshifted sources as compared to the redshifted ones. 

3. Non-Copernican Models 

We now examine non-Copernican models where the observer is close to the center of a 
very large underdensity, which have been shown to mimic accelerated expansion without 
dark energy* We will restrict our discussion to void models based on the LTB metric 
and we will consider both compensated and uncompensated profiles. The former are 

* In [71, 72] it has been proposed that a model with a central overdensity can reproduce the lightcone 
expansion rate and matter density of the concordance ACDM model. However this is possible only by 
invoking a strongly non-simultaneous big bang, and it is not clear if this can be made compatible with 
other observations, such as the CMB, as we briefly discuss later in Section 3.1. 
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very similar to the Swiss-cheese holes of the previous Section, the only difference being 
the much larger size and the central position of the observer. Uncompensated voids, 
instead, are solutions that approach a background FLRW model only asymptotically. 
Our discussion will be qualitatively valid for void models not based on the LTB 
metric, such as the so-called "Hubble-Bubble" models, in which the inner underdensity 
and the outer homogenous region are simply described with two disconnected FLRW 
metrics [73, 74, 75]. 

3.1. SNe Observations 

In the past 15 years it has been extensively studied (see, for example, [17, 18, 76, 
77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 71, 90, 91, 92, 93, 94, 95]) how an 
observer inside a spherical underdensity expanding faster than the background sees 
apparent acceleration. This effect is easy to understand: standard candles are confined 
to the light cone and hence temporal changes can be replaced by spatial changes along 
photon geodesies. In this case, "faster expansion now than before" is simply replaced 
by "faster expansion here than there" . This is why a void model can mimic the effect of 
a cosmological constant if it extends to the point in space/time where the dark energy 
becomes subdominant: a typical scenario that can mimic the late-time acceleration of 
the concordance ACDM model consists of a deep void extending for 1-3 Gpc. 

We may understand what happens by looking again at Eq. (1). In a void model 
slightly smaller than the horizon the dominant term in the Newtonian gauge is the 
Doppler effect, which induces a large net correction to the redshift and modifies the 
di(z) curve in the same way of an accelerating universe; in other words, it is as if all 
sources in the underdense region had a radial collective peculiar velocity which gives 
them an extra redshift. This is the dominant effect at low redshift and therefore the 
most relevant for SNe. This behaviour has been extensively studied through analytical 
approximations and both in the synchronous and Newtonian gauge in [46, 47, 83], an 
example of which is given in the left panel of Fig. 5. In addition there is a correction 
to the distance itself, which, however, becomes important only at high z: a change in 
the angular metric element Y(r, t) induces lensing corrections to the area under which 
a single source is seen. 

More formally one may consider the freedom available within LTB models: as 
explained in Appendix A, the LTB metric features indeed two free functions. In the 
usual description the LTB model is analyzed in comoving and synchronous coordinates, 
where one free function called the "bang time" function te(r) sets the time of the big 
bang at each value of the coordinate r, and the other function can be taken as the 
density profile or the r-dependent curvature (depending on arbitrary redefinitions of the 
radial coordinate, which leave the LTB metric in the same form). It is enough to adjust 
one free function (e.g., the curvature) to obtain any distance-redshift relation, (j as stated 

jj Strictly speaking it is not possible to have a negative deceleration parameter go a t the coordinate 
point of the observer, unless one uses a density profile which is spiked at the origin [96, 97]. Nonetheless, 
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Figure 5. Left panel: comparison between analytic (short-dashed line) and numerical 
(solid line) cLl{z) curves as obtained in [83]. Also plotted are the EdS curve (long- 
dashed line) and the ACDM curve with J7a = 0.7 (dotted line). Right panel: SN 
data (SDSS SN) and predictions normalized to a reference model with VLm = 0.3 and 
n K = 0.7 (OCDM) for the best-fit model of [92] (CMB + BAO + SNc + HST). Also 
plotted are the curves relative to ACDM model (dashed green) , void model with curved 
background (dotted blue) and void model with flat background (fine dotted magenta) . 



by Ref. [77] and illustrated, for example, by Refs. [85, 89] where the luminosity distance 
of the concordance model is reproduced (see also [98]). By adjusting the other free 
function, it is possible to obtain also the light-cone matter density (or galaxy number 
count) of the concordance model [89]. However, it is probably preferable to avoid an 
inhomogeneous bang function (see the simultaneous big bang condition in Eq. (A.4)), 
which would introduce very large inhomogeneities in the past, strongly at odds with the 
inflationary paradigm [46, 99]. It is possible, however, to demand that t' B (r) ~ 0, but not 
strictly zero. In this case one should guarantee that observations are compatible with 
the inhomogeneities in the big bang time. In particular, constraints can be imposed by 
demanding that, for a radius Tlss and time tiss which correspond to the Last Scattering 
Surface, the density perturbations associated to t' B (r) are small, in order not to spoil 
CMB observations [100]. This implies a strong suppression at present time on t' B at 
r ~ r Lss, since these are decaying modes. More freedom could be allowed for r <C r^ss, 
keeping in mind however that one should not introduce too large secondary effects on 
the CMB, and not spoil galaxy number counts and galaxy ages. In the rest of the review 
we will always set the bang function to ts = 0, which is the simplest choice and also 
makes the model more constrained with respect to the t' B ~ case. 

We show in the right panel of Fig. 5 an example from [92] of how SNe data can be 
fitted with a void model. The top left panel of Fig. 6 illustrates the fact that concordance 
and void model equally well fit the SNe distance-redshift diagram, showing confidence 
level contours on the void radius ro and f^A.out for the ALTB model of Ref. [95] (the 
subscript "out" labels quantities relative to the background model). In [95] the void 
has the smooth density profile already used in [84] , which depends crucially only on the 

it is easy to evade this constraint observationally with a rapidly varying function or by smoothing the 
density profile on a appropriately small scale [85], since it is true only in one point. This fact, however, 
might be of some observational interest for future precise measurements at very low redshifts. 
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Figure 6. 1, 2 and 3a confidence level contours for the asymptotically flat ALTB 
model of [95] on the void radius tq and f^A,out 

with n s = 0.96 and t Q = 13.7 Gyr. The 
four smaller panels on the left show the contours for the independent likelihoods per 
observable, while the larger panel on the right shows the contours for the combined 
likelihood. In the panel representing measurements of the local Hubble constant, the 
results relative to i?so6 °f Eq- (8) are shown as filled contours, while the ones relative 
to -£/ro9 °f Eq. (9) are shown as lines. The same labelling holds for the panel relative 
to the combined observables. In the panel relative to CMB constraints we also show 
confidence levels for the likelihood marginalized over n s (dot-dashed contours). In the 
panels, the x-axis represents pure-matter void models, while the y-axis the standard 
FLRW model. 



void radius and depth, the two main physical quantities describing an underdensity. 
3.2. CMB 

Dust LTB void models are matched to the background metric at a redshift at which 
radiation is still negligible; a value of z ~ 100 satisfies, for example, this requirement. 
In this way the Last Scattering Surface (LSS), which is responsible for most of the CMB 
anisotropies, is outside the inhomogeneous patch (and so there are no Doppler effects) 
and a standard analysis of the primordial CMB power spectrum is possible (but perhaps 
not of the BAO as we will discuss in the next Section). Moreover, the lensing effects 
discussed in the previous Section are almost vanishing, as light rays always go radially 
for a central observer in spherically symmetric void modelsf f and the angular diameter 
distance <Ia{z) = Y(r(z),t(z)) is very close to the one of the background model: the 
scale factor Y exactly matches the background one (Y = ar for r > r ) and the only 
difference is in a small correction (redshift effect) to the arguments r(z) and t(z) when 
solving for the geodesies. This small shift in the geodesies slightly changes the distance 



ff Light propagation is qualitatively different for the non-sphcrically symmetric setup of Fig. 1. 
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to the LSS and it can be reabsorbed into a redefinition of the photon temperature of the 
CMB. In other words, everything can be analyzed with an effective FLRW model with 
a slightly different value of To [86, 92], which will generally differ from the usual 2.726 
K and consequently, the LSS is located at a slightly different redshift. Note also that 
there are other contributions to the CMB coming from secondary effects, which may 
differ from FLRW and are due to the photons traveling through the inhomogeneities 
inside the void, such as the ISW effect or lensing; in the literature these effects have 
not been considered, since they are subdominant and since studying them would require 
knowledge of the growth of perturbations in an LTB model, which is not yet completely 
understood (see, however, [99, 101]). 

Now, in a compensated model the correction to T is of order (Zo/^hor) 2 [92], therefore 
small for void radiuses of about 1-3 Gpc (0(1%), [86, 92, 95] ), and since it is well-known 
that it is not possible to fit the CMB observations with an EdS model, it follows that 
(assuming the standard primordial power spectrum) compensated void models which 
are asymptotically flat are ruled out. This holds independently of how the LTB free 
functions are adjusted and is shown by the bottom left panel of Fig. 6 taken from [95]. 
Similarly, the asymptotically flat void model has been shown by [92] to be excluded 
with a Ax 2 ~ 35. In [95] the CMB power spectrum has been analysed using accurate 
fits for the positions of the first three peaks and the first trough and the relative heights 
of the second and third peak relative to the first one, while in [92] the full C\ spectrum 
has been studied using a modified version of C0SM0MC. The agreement between the two 
approaches strengthens their conclusions. 

One way of improving the fits is to consider uncompensated voids, which have a 
large correction to the monopole temperature T and therefore a different distance to 
the LSS. It has been shown [86, 92] in fact that such a model with an (9(20% — 40%) 
correction to the monopole, taken at face value, can have a much better fit than the 
asymptotically flat model. It is not clear, however, if in this case it is consistent to 
analyze the CMB in the standard way, since the LSS itself is affected by a radial peculiar 
velocity, and so these models may have to be considered only illustrative of the potential 
effects they cause to the observables [92]. 

A simpler way out of the shortcomings of the asymptotically flat model is to consider 
void models that are asymptotically curved [92, 95] and so also with a different distance 
to the LSS: a closed model (e.g., Qk ~ —0.2) older than the ACDM concordance model 
(e.g., to ~ 16 -T- 18 Gyr) can indeed fit the observed CMB, as shown by Table 1 which 
corresponds to the profiles of Fig. 7 studied in Ref. [92]. Note that these fits include also 
the value of H and BAO constraints, as discussed in the next two Sections. We can 
summarize the previous discussion by stating that CMB observations fix the background 
metric of a compensated void model, with a weak dependence on the LTB free functions. 
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Modpl 


CMB 


BAO 


SNe 


^soe 


Total v 2 


ACDM 


3372.1 


3.2 


239.3 


0.4 


3615.0 


Profile A (Curved Void) 


3377.1 


4.3 


240.7 


6.6 


3628.7 


Profile B 


3377.2 


0.6 


235.3 


5.1 


3618.2 


Profile C 


3376.9 


1.0 


234.9 


3.7 


3616.5 


Profile D 


3376.7 


3.8 


233.9 


2.2 


3616.6 


Profile E 


3372.9 


3.2 


241.9 


1.1 


3619.1 



Table 1. A breakdown of the total \ 2 for each dataset, for fitting simultaneously to 
CMB + BAO + SNe + H sm . The profiles refer to Fig. 7 and were studied in Ref. [92]. 



3.3. Local H 

It is well known that a curved FLRW model can fit the CMB only at the price of a very 
low value of the Hubble rate and this remains true, as discussed in the previous Section, 
for the Hubble parameter H out of the FLRW region outside the void. When comparing 
with observations, however, the presence of the void itself alleviates the problem as the 
observer experiences a higher local expansion rate Hq. The jump AH = 100 Ah km 
s _1 Mpc -1 between H out and H Q is constrained by a good fit to the SNe to a value of 
Ah ~ 0.2 and so we can conclude that CMB observations fix the magnitude of the local 
Hubble rate, which can then be confronted with the observations. We will consider the 
following two results from Ref. [102] and Ref. [103]: 

Hsoe = 62.3 ± 5.2 km s^Mpc -1 Sandage et al. 2006 , (8) 
H Rog = 74.2 ± 3.6 km s^Mpc -1 Riess et al. 2009 . (9) 

The general trend is that the CMB favors values of the local expansion rate which are 
at most of about Hq w 50 km s _1 Mpc _1 [92, 95], so in disagreement with i^Rog, but in 
marginal agreement with Hsoe- It is, however, possible to exploit further the freedom 
of the LTB models in order to increase the local value of H by adjusting the LTB 
free functions in order to have a small local patch around the observer, with a higher 
expansion rate, that does not disrupt the SNe fit [92] and this has been shown to also 
give a better overall combined fit to cosmological observables, as can be seen in Table 1. 

Another possibility considered in the literature, alternative to the simplest 
asymptotically flat models, is to adopt a non-standard (but physically motivated) 
primordial power spectrum as suggested by Ref. [104] (see also [94]), where a void model 
matched to EdS was shown to be in agreement with SNe, CMB and H observations. Yet 
another possibility is to consider an inhomogeneous profile for the radiation [93], thus 
allowing for the extra freedom necessary to accommodate the CMB power spectrum. 

3.4- BAO and Power Spectrum 

A baryon acoustic peak has been detected from SDSS and 2dFGRS galaxy catalogues at 
a mean redshift of 0.2 and 0.35 [105]. This measurement constrains theoretical models 




Figure 7. Curvature (top) and density (bottom) void profiles studied by Ref. [92], 
whose corresponding x 2 are displayed in Table 1. By modeling the void with the 
profiles B and C it is possible to obtain an excellent fit to the BAO observations. 



by means of the following ratio: 

Ls 



(10) 



dy{z) 

where Ls is the comoving sound horizon scale at recombination and dy is a combination 
of angular and radial distance defined as follows: 

d v = [(l + z) 2 d 2 A d z } 1 / 3 . (11) 

Void models able to fit the SNe have to extend at least up to z ~ 1 and so the BAO 
feature is well inside the inhomogeneous patch, where spatial gradients are strong. It is 
possible to compute in a given LTB model [92] the above distances, taking into account 
that the transverse and the radial expansion rates are different. However, in order to 
compare with observations one needs to make assumptions, which make the validity of 
the analysis unclear. The first assumption is to ignore the effects of the background on 
the growth of perturbations, which is unlikely to be correct in an LTB model because 
of the anisotropics in the expansion rate. As for the secondary effects on the CMB 
discussed in Section 3.2, this has not yet been addressed in the literature since the 
growth of perturbations in an LTB model is not clearly understood (see Ref. [99, 101]). 
The second assumption is to consider that Ls is given by FLRW recombination physics, 
which is unlikely to be correct in this kind of models: it is true that the baryon acoustic 
scale is imprinted in the sky at early times, when the universe was supposed to be more 
homogenous, but it has been imprinted in a spatial location in which the void itself 
was going to develop and therefore this should be analyzed in a radiation dominated 
era within a spherically symmetric background. Further, in order to perform a more 
accurate analysis one should compute 9(z) for a generic z (and not only for the two 
values above) and then weight it, for instance, by the corresponding number density of 
observed objects, especially if the void profile is rapidly varying in this redshift range [92]. 
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If one insists, with the above caveats, in confronting LTB models against BAO 
observations, one finds that a flat (near the origin) density or curvature profile generally 
has difficulties in fitting the BAO scale, as shown in the example of Fig. 6. It is however 
possible by a more specified tuning of the void profile to fit the BAO data without 
cosmological constant. A better fit than the ACDM model can be found indeed by 
adopting the profiles shown in Fig. 7, as can be seen in Table 1 from [92]. 

Similar considerations apply to other constraints coming from the matter power 
spectrum. Here also it is essential to study the growth of perturbations in an LTB model, 
and again this has not yet been done (see however [106] for recent improvements using 
numerical simulations). A fit of these data can nonetheless be performed under some 
rough assumptions, in order to get some qualitative indication. In [92], for example, fits 
have been performed assuming the growth of structures as in an approximate effective 
FLRW model built up from the parameters at the centre of the void. The results show 
that in this case the void gives a bad fit, be it with flxput = or with f^ iOUt ^ 0, since 
matter power spectra favour a value for the combination fioM.in^ of about 0.2, while 
the curved void from the fit with CMB + BAO + SNe + -£^so6 has a best fit value of 
about 0.09. Note, however, that the fit is significantly better for the LRG subset [107] 
than for the SDSS main sample [108]. In fact, the void does predict too much power on 
large scales and the LRG data prefer this, as compared to the SDSS data. 

3.5. Constraints on the observer's position. 

Most works in the literature take the observer to be at the center of the LTB void in 
order to recover the observed isotropy of the universe. This is undoubtedly a heavily 
fine-tuned configuration, especially for the very large voids considered. It is therefore 
interesting to see which actual limits are placed on the off-center position of the observer 
by present or future observational data. 

Using the anisotropies in the luminosity of the observed supernovae (see [109] for a 
claim of detection of discrepancy between the equatorial North and South hemispheres), 
it has been found [110, 111] that SNe observations constrain the observer's position to 
about 200 Mpc from the void center. Much tighter constraints come, however, from the 
observed dipole of the CMB; this is best understood in the Newtonian picture, see [112] 
for exact numerical results and [113] for analytical results valid for general spherically- 
symmetric spacetimes. An observer close to the center and comoving with the LTB 
metric has a "background" peculiar velocity with respect to the background model 
(and the void center) of Av = AH <i obs , which is completely due to the inhomogeneous 
expansion rate (as opposed to the "real" peculiar velocities in the standard FLRW 
paradigm; see Ref. [30] for a discussion about peculiar velocities in inhomogeneous 
universes). This causes a Doppler shift in the CMB temperature and a consequent 
dipole of magnitude [112]: 




(12) 
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As said earlier, SNe data typically constrain void models to have Ah ~ 0.2 and so the 
dipole of Eq. (12) has a magnitude of the order of the observed one for an observer 
displaced from the center of about 20-30 Mpc. For smaller displacements, the LTB- 
induced dipole is smaller and peculiar velocities are invoked, as in the FLRW model, to 
match the observed value. For larger displacements, the induced dipole is too big and 
the observer's background peculiar velocity has to be compensated by a real peculiar 
velocity directed towards the center of the void. It has been found by [114, 115] that, for 
typical values of real peculiar velocities, the position of the observer is constrained to be 
up to 60-80 Mpc from the center, which is about few percents of the typical radius of 
a void model. The fine tuning on the observer's position is therefore heavy. It is worth 
mentioning, however, that such a tuned position of the observer could simultaneously 
explain [92] the recent observation of large bulk flow velocities [116] of order 600 -j- 1000 
km/sec in a local region of about 300/i _1 Mpc close to the direction of the observed 
CMB dipole. An observer displaced from the center would indeed observe a difference 
in velocities in the sky with a dipolar shape aligned with our local CMB dipole, with a 
magnitude of the order of Av in agreement with [116]. 

Finally, it has been shown by [114, 117] that tighter constraints will come from 
future high-precision astrometric observations of distant quasars. Cosmic parallax 
measurements can indeed put stringent bounds (~10 Mpc) on the off-center position of 
the observer, which are independent from the constraints relative to the CMB dipole. 
In particular, the effect depends on the source and observer peculiar velocities, and this 
may break the degeneracy between real and background peculiar velocities discussed 
above. Cosmic parallax together with the redshift drift (to be discussed in the next 
Section) belong to the realm of the so-called "real-time" cosmology [118]. 

It is important to stress that the previous discussion is relative to the LTB metric in 
which only one center position can be meaningfully defined. All the previous constraints 
could indeed be in principle alleviated if one considers non spherically-symmetric metrics 
which allow for more (or none) center positions [42]. In other words, in more general 
metrics the center location of the LTB models could become an extended volume, thus 
reducing the fine tuning required on the observer's position. It would be interesting 
to quantitatively define the amount of fine tuning within a void model by taking, for 
example, the ratio of the volume occupied by observers who see a CMB dipole smaller 
than or equal to the observed one to the total volume of the void region. In this way 
one could precisely see how weakened is the fine tuning on the observer's position in 
more general geometries as compared to the LTB case. 

3.6. General tests of the Copernican Principle 

We will now discuss more general tests and observables able to further constrain the 
void models, which apply also when the observer is exactly at the center. 

Even if we do live at the center of the LTB void, other halos and clusters in the 
universe do not and the corresponding observers will see a large dipole in the CMB. Such 
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a dipole would manifest to us observationally through the kinematic Sunyaev-Zel'Dovich 
(kSZ) effect. The hot electrons inside a cluster distort indeed the CMB spectrum 
through inverse Compton scattering, in which the low energy CMB photons receive 
energy boosts during collisions with the high-energy cluster electrons. In the (first-order) 
thermal Sunyaev-Zel'Dovich effect the CMB photons interact with electrons that have 
high energies due to their temperature, while in the (second-order) kinematic Sunyaev- 
Zel'Dovich effect the CMB photons interact with electrons that have high energies due 
to their bulk velocity. The kSZ effect is proportional to the radial velocity of the cluster 
with respect to the CMB scattering surface and can be employed to map the cosmic 
peculiar velocity field inside our own light cone. 

Ref. [119] by examining available kSZ measurements was able to put bounds on 
the size of the void radius, which is then constrained to be smaller than about 1.5 
Gpc. They also conclude that, within their adopted void modeling, the kSZ constraints 
cannot be satisfied if the simultaneous big-bang conditions is demanded. This clearly 
shows how tests beyond the observer's lightcone are crucial to constrain void models (see 
also [120]). Moreover, it would be interesting to quantify the amount of inhomogeneity 
in the bang function (and so of decaying modes) that is necessary to fit the kSZ data. 
Ref. [121] examines the same dataset considered by Ref. [119] and finds their void model 
ruled out due to the high peculiar velocities. The different conclusion with respect 
to Ref. [119] might be attributed to the different modeling of the void; see indeed 
Fig. 3 of [119] and Fig. 4 of [121], which show that the two models have very different 
radial velocities at z > 0.5. Ref. [121] considers also the case of an inhomogeneous 
decoupling hypersurface. A large kSZ effect is indeed due to a large velocity between 
the comoving clusters and the CMB frame and this can be alleviated if radial non- 
adiabatic (isocurvature) inhomogeneities in the non-relativistic matter on the decoupling 
hypersurface are introduced. More work is necessary to understand the implications of 
such a non-standard universe (see also [93]). 

Strong results have been obtained recently by Ref. [122] where the kSZ effect due 
to all free electrons was calculated and confronted with the limits coming from the 
measurements of the small scale anisotropy power spectrum by the Atacama cosmology 
telescope. According to Ref. [122] void models are robustly excluded. This test is 
worth careful consideration, in order to conclude on the possible definitive failure of 
the void models, examining in detail its assumptions and calculations (see Ref. [123] of 
this Special Issue). For instance, one caveat is that it is assumed that the growth of 
perturbations is well described by a ACDM model and the validity of this assumption is 
unclear. Moreover, non-adiabatic models could also help in avoiding these constraints. 

Other tests of void models (and so of the Copernican Principle) have been 
considered in the literature [124, 125, 126]. Rescattering of photons by off-center 
reionized structures can distort, for example, the CMB blackbody spectrum via 
Compton y-distortion [75, 94]. 

Another potentially interesting observable is the redshift drift, namely the temporal 
variation of the redshift of distant sources like quasars as a tracer of the background 
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cosmological expansion, which was first discussed by Sandage in 1962 [127] (see also 
[128, 129]). It relies on high-precision spectroscopy and the necessary statistical 
sensitivity could be reached with the next generation of optical telescopes. In particular, 
a redshift drift measurement A t z over a time-span At could distinguish between the real 
acceleration driven by dark energy (A t z > 0) and the apparent acceleration of the void 
models (A t z < 0). Most importantly, it is a direct measurement of the local expansion 
rate of the universe which is independent from the evidence for acceleration given by the 
SNe, thus being an important test independent of the calibration of standard candles 
and the relative uncertainties. The redshift drift in relation to LTB void models has 
been examined, for example, in [90, 114, 130, 131]. 

Finally, one can further constrain the LTB void models by means of other 
observational quantities such as source number counts [98] and, more recently, galaxy 
ages [132]. 

4. Discussion and outlook 

In this paper we have focused on the modeling of voids, which are the characterizing 
feature of the late inhomogeneous universe. We have in particular discussed the 
possibility that the observed dark energy could be explained away by the effect of 
large-scale nonlinear inhomogeneities. We considered both the Copernican and non- 
Copernican paradigm and we discussed how cosmological observations can constrain 
these models. Work still has to be done in order to thoroughly explore inhomogeneous 
models, but it is nonetheless worth drawing some partial conclusions and maybe indicate 
possible future directions. 

Within Copernican models it is useful to distinguish between two non- exclusive 
effects of large-scale inhomogeneities on cosmological observables. The first one 
addresses the question of how the cosmological background reacts to the nonlinear 
structure formation, while the second focuses on the propagation of photons in the 
inhomogeneous universe. The former effect has been called sometimes "strong" 
backreaction and the latter "weak" backreaction (see Ref. [30] for the precise definitions). 
Strong backreaction has been discussed in other contributions to this Special Issue [22, 
23, 24, 25, 26, 27, 28], but we can nonetheless say that a common agreement on its 
magnitude has not yet been reached. The models reviewed in this paper cannot answer 
this question because by construction the evolution of the background is unaffected by 
the inhomogeneities introduced, despite them being fully nonlinear. On the other hand 
we have studied photon propagation in these models and, from what we have presented 
here, it seems that a single (weak backreaction) effect on photons is not sufficient to 
have a paradigm shift in which dark energy is no longer needed [38, 39, 46, 47]. We 
think, however, that Copernican models should be studied in more detail, because even 
if each single effect is small, of order 1%-10%, their sum might not be, especially taking 
into account selection effects. This conjecture is strengthened by the fact that different 
effects of inhomogeneities seem to pull into the same direction with less and less need 
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for dark energy [61]: to conclude on their viability as the possible explanation of the 
observed acceleration, it is therefore crucial to consider all possible effects together in 
as realistic a model as possible. 

With non-Copernican models the situation is, to some extent, opposite. Since the 
very first works of 10-15 years ago it was indeed possible to mimic the dark energy, and 
the issue of the viability of these models was more philosophical than observational. 
However, since then, all the research has focused on exploiting the freedom available in 
order to accommodate more observables and recent data. Most literature has focused 
on models which employ only one of the two free functions in LTB, setting the other by 
asking a homogenous bang function (relaxing this assumption to an almost simultaneous 
big bang could make it easier to satisfy the constraints). With more data and more 
research it appears that for one single free function, it is by now not sufficient to consider 
the simplest void models embedded in EdS to fit all data, but it is necessary to introduce 
some features, like overall background curvature, more elaborated profiles, or nontrivial 
primordial spectra for the CMB. In this way, at least when considering CMB+SNe+ifo; 
such models succeed in fitting the data with more or less success, depending on the 
specific profiles [92, 95]. The inclusion of the BAO observations can be done only at 
a qualitative level, as it would require knowledge of evolution of perturbations in LTB 
during the radiation and matter eras; the present analyses have done under simplifying 
assumptions and indicate that these data could be also fit by LTB, especially with 
peculiar shapes of the density profile. Similarly, one may also attempt a fit of the 
matter power spectrum under the same assumptions, while the result in this case is 
that the fit is significantly worse [92]. 

Finally, a significant and important amount of work has shown that these models 
can be severely constrained using the fact that there are large radial peculiar velocities. 
The position of the observer is indeed constrained to be at a distance of at most about 
60-80 Mpc from the center of the void, in order to avoid a too large dipole in the CMB. 
This represents a strong fine-tuning on the observer's position (which perhaps could be 
weakened by modeling the void with a less-symmetric metric than LTB, as we briefly 
discuss in 3.5), but on the other hand an off-center observer could explain the recently 
claimed large bulk flows on very large scales [116]. A second very interesting set of 
constraints comes from the fact that light from CMB is rescattered by electrons which 
live in structures with a large radial collective velocity. These tests, which apply also 
if the observer is exactly at the centre, put severe constraints on such models and may 
turn out to rule out most, if not all, of the void scenarios. Some of these constraints 
rely again on simplifying assumptions on the growth of perturbations, and therefore 
more refined analyses might be needed to conclude on the viability of non-Copernican 
models. 
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Appendix A. LTB basic formalism 

We will now quickly review the conventional LTB formalism. For more details see, for 
example, Ref. [46, 47, 81, 40]. The line element of the spherically symmetric LTB model 
in comoving and synchronous coordinates can be written as (c = 1): 

ds 2 = -dt 2 + ^T'^ dr 2 + Y 2 (r, t){dd 2 + sin 2 Odcj?) , (A.l) 
1 + lEyr) 

where Y(r,t) is the scale function, the prime denotes derivation with respect to the 
coordinate radius r and the arbitrary function E(r) is the curvature function which is 
related to the spatial Ricci scalar by 1Z = —A(EY)'/(Y 2 Y') [51]. Throughout the paper 
we will loosely refer to E when talking of spatial curvature in the LTB model. The 
FLRW solution is recovered by setting Y(r, t) — > a(t) r and E(r) — > —k r 2 /2 throughout 
the equations. Note that in the LTB space the transverse expansion rate = Y/Y 
will generally differ from the longitudinal expansion rate Hl = Y'/Y'. 

The dynamics of the model is governed by the following equation [34]: 

Y 2 2F(r) SttG 2E(r) 

Y2 = ~Y^~ + ~1~ Pa + ~Y*~ ' ( > 

where the dot denotes derivation with respect to the coordinate time t and p\ = A/8nG 
is the energy density associated with the cosmological constant. The arbitrary function 
F(r) (actually a constant of integration) represents the effective gravitating mass and 
is related to the local dust energy density pM(r,t) through the equation AnG pm = 
F'/(Y 2 Y r ). It is interesting to note that in the last equation the curvature term is 
missing and so the gravitating mass - which is the quantity that matters in matching 
the metric to the background - differs from the invariant mass [34]. This leads to 
"strong" backreaction [30] effects in LTB models as it implies that the expansion rate of 
a given shell (e.g., the border of the hole) is not sourced by the averaged mass within. 
We will discuss this in Appendix B for the case of Swiss-cheese models. 
It is useful to rewrite Eq. (A. 2) in the following, more familiar form 

where H (r) = H T (r,t ), Y (r) = Y(r,t ) and the (present-day) density parameters are 



H 2 (r,t) = H 2 (r) 



, (A.3) 



2F(r) _ 8ttG Pa 



H 2 (r)Y*(r) ' ^ ' ~ 3 H { 



r 
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2E(r) 

Eq. (A. 3) can be used to determine the age of the universe at a radial coordinate r: 

1 

to - t B (r) = [ dX . (A.4) 

H (r) J ^tt M (r)x- l + tt K (r) + tt K (r)x 2 



o 



One can constrain the models by requiring a simultaneous big bang, i.e., by setting 
^bO") — 0. Simultaneous big bang excludes decaying modes which would be strongly 
in contradiction with the inflationary paradigm [46, 99]. Note, however, that an almost 
simultaneous big bang t' B {r) w could be used [71] as discussed in Section 3.1. 

LTB models feature three arbitrary functions. Within the present formalism they 
are taken as Q M (r), t B {r) and Y (r). One of these is but an expression of the gauge 
freedom, which we may fix by setting Yq{t) oc r or, equivalently, F(r) oc r 3 . 

Appendix B. "Strong" backreaction in exact Swiss-cheese models 

We will now argue that "strong" backreaction [30] effects on the background evolution 
are small for realistic Swiss-cheese models based on LTB metrics which feature a void 
at the center. See Ref. [40] for a discussion of more general LTB setups. 

Let us start by calculating the average expansion rate of an LTB hole at a fixed time 
t: any deviation from the background value H will signal strong backreaction effects 
due to the matter inhomogeneities. This is formally done by averaging the expansion 
scalar 9 = Hi + 2Ht- It is easy to see that: 



= f v 6dV = fr°dr6Y 2 Y'/VTT2E = Ve%i Y 
W V ) r Q ° drY 2 Y>/^/T+2E V Y 



3H, (B.l) 



where V is the proper 3-dimensional volume of the space slices, and we remind that ro 
and Iq are comoving and proper radiuses of the hole, which is taken to be much smaller 
than the horizon: Iq <C l^ OT . Eq. (B.l) shows that if the curvature E is small then 
the average expansion rate (6)/ 3 is close to the background value H and the strong 
backreaction is negligible. This is indeed the case of flat LTB models (E(r) = 0) for 
which it is easy to show that the kinematical backreaction Qd of Ref. [21] is exactly 
zero no matter how the free functions are specified [51]. 

We have now to estimate the order of magnitude of E for realistic LTB models with 
a void at the center. Let us consider for simplicity an EdS background and fix the gauge 
by setting Yq(j-) oc r. It is useful to rewrite Eq. (A. 2) as: 

^(r) = ^(r,t)r 2 (r,t)-^|, (B.2) 

which can be understood as that the total energy per unit of mass of the shell r is given 
by its kinetic energy per unit of mass (iff. Y 2 = Y 2 ) plus its potential energy per unit of 
mass due to the total gravitating mass up to the shell r. Note that thanks to spherical 
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symmetry one is able to define a potential energy also in cases far away from nearly 
Newtonian ones and that the potential energy is related to the curvature [34]. 

The curvature function E vanishes at the center (generally it is E oc r 2 [77]) and at 
the border where the hole is matched to the EdS metric, for which E = and so the two 
terms of Eq. (B.2) cancel each other: see, for example, Fig. 4 of Ref. [38]. Moreover, at 
the center of the hole there is a void and so the second term in Eq. (B.2) is in magnitude 
< as compared to the first one and so we can evaluate E by considering the first term 
only. Roughly, for realistic voids, it is Ht ~ H and so we have: 



that is, the curvature is indeed small for sub-horizon holes [30] (see also [133]). In 
Eq. (B.3) we have evaluated Y by taking its maximum value Y(ro,t) = Zo: we remind 
indeed that Y = at the origin and Y' > as Y' = would give shell crossing (we 
are assuming F'(r) > 0). We stress that even if the compensating overdense shell 
features matter contrasts ^> 1, still F(r) is close to the EdS value as it is an integrated 
quantity which smoothly tends to the background value. This shows how spherical 
symmetry severely restricts the effective inhomogeneities allowed in an LTB model. As 
said earlier, the dynamics of the shell r depends only on the total mass within and does 
not depend on the mass outside, that is, density inhomogeneities are already averaged 
out as far as the dynamics of the shell is concerned (similar considerations apply to the 
curvature function E which also is an average quantity as far as the spatial Ricci scalar 
is concerned). 

Eq. (B.3) connects the curvature to the expansion rate in the void. It is possible to 
evade the latter constraint by considering LTB models which feature Ht 3> H. These 
models, however, can only be matched to the background metric asymptotically and not 
at a finite radius Z as with Swiss-cheese models. Such a strong expansion rate, besides 
being unrealistic, would indeed cause shell crossing in the overdense shell surrounding 
the void in a very short time interval. As we have said before, the dynamics of the inner 
void shells is unaffected by the outer dense shells which will be (no matter how the free 
functions are adjusted) squeezed towards the border of the hole causing shell crossing 
in a time At ~ (l - l /2)/(AH l /2) ~ AH' 1 ~ H^ 1 < H~ l , where we estimated 
the edge of the void to be at Z /2 (and so its peculiar velocity to be AH Z /2) and the 
"space" to cover to be Z — Zo/2. See Ref. [134] for a discussion of the evolution of general 
LTB profiles. 

We have shown, at least qualitatively, what causes the curvature to be small: 
AH ~ 1 on sub-horizon scales Zo- in other words, the curvature is small if so is 
the quantity AH Z , which in Ref. [30] was interpreted as a background peculiar 
velocity. As said before, it is $ oc E [47] and one could re-interpret the smallness 
of the strong backreaction as the fact that the inhomogeneous universe (in this case the 
Swiss-cheese model) can be described by means of a newtonian perturbed metric with 
a small potential $. The calculations of this Appendix, however, show that this is an 
assumption of the model and not a general feature of the inhomogeneous universe. It 




(B.3) 
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is indeed clear that in the present case the curvature (and so the newtonian potential) 
is small because the (sub-horizon) inhomogeneities are matched to an a priori chosen 
background: in other words the strong backreaction is assumed small since the beginning 
by demanding small background peculiar velocities [30, 135]. 
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